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H ' We report on a search for Ha absorption in four exoplanets. Strong features 

at Ha are observed in the transmission spectra of both HD 189733b and HD 
209458b. We attempt to characterize and remove the effects of stellar variability 
pq ' in HD 189733b, and along with an empirical Monte Carlo test the results im- 

ply a statistically significant transit-dependent feature of (—8.72 ± 1.48) x 10~^ 
Ph ■ integrated over a 16 A band relative to the adjacent continuum. We interpret 

. ' this as the first detection of this line in an exoplanetary atmosphere. A previous 

p ■ detection of Lya in HD 189733b's atmosphere allows us to calculate an exci- 

Q ■ tation temperature for hydrogen, T^xc = 2.6 x 10^ K. This calculation depends 

C^ ■ and argue that T^xc is very likely much greater than the radiative equilibrium 



> 



temperature (the temperature the planet is assumed to be at based on stellar 
radiation and the planetary distance) of HD 189733b. A large Texc implies a 



Tj- ■ very low density that is not in thermodynamic equilibrium the planet's lower 

atmosphere. We argue that the n = 2 hydrogen required to cause Ha absorption 
"^ ■ in the atmosphere is created as a result of the greater UV flux at HD 189733b, 

^ ■ which has the smallest orbit and most chromospherically active central star in 

^ ■ our sample. Though the overall integration of HD 209458b's transmission spec- 

trum over a wide band is consistent with zero, it contains a dramatic, statistically 
significant feature in the transmission spectrum with reflectional symmetry. We 
rS ■ discuss possible physical processes that could cause this feature. Our remaining 

c^ . two targets (HD 147506b and HD 149026b) do not show any clear features, so 

we place upper limits on their Ha absorption levels. 
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INTRODUCTION AND BACKGROUND 



Extrasolar planets that are known to transit their central star present an important op- 



portunity for astr onomers to study their atmospheres. Fo 



siting exoplanet fjCharbonneau et al.l l200d : iHenrv et al 



lowing; the first detection of a tran- 



2000 



Mazeh et al.ll2000f ). the first 



atmosphere of an exoplanet was discovered by lCharbonneau et al.l (120021 ) . who found sodium 



absorption in the atmosphere of HD 209458b. In the en suing years there have been addi 



tiona l detections in exop lanetary at mospheres of sodium (iRedfield et al.ll2008l : [Snellen et al. 



20081 : IWood et al.ll201ll ). potassiur n (IColon et a 



cludi r ig H^O, CH^, CO, and C O9 (IBarman 



2008; Swain et al. 



2008, 2009 



m 



I2010 



2007; 



2010 



Beaulieu et al. 



Tinetti et al. 



2008 



Sing et al.l 120111). and molecules in 



2010; Grillmair et al. 



20071)111 These detections are presumed 



to come from the bound atmosphere of the host planets. 

All atmospheric detections can be thought of as wavelength-dependent variations in 
apparent planetary radius, observed using either narrow-band photometry or spectroscopy. 
With information about the absorbing species, this can be translated into a height for the 
absorbing material, which may indicate that the absorbing material extends to the unbound 



portion of an exoplanet's atmosphere, i.e., the exosphere. T 
exoplanets, with HD 209458b being the most w ell-studied. 



lis has been observed in multiple 



Vidal-Madiar et all f l2003l . l2004i ) 



found H I, C II, and O I in HD 209458b, while JLinsky et al.l (J2010l ) confirmed the presence 
of H I and detected Si III. ISchlawin et al.l (120 lOh al so observed Si IV in this target. O ther 
exospheric detections include H I in HD 189733b bv iLecavelier Pes Etangs et al.l (120101 ) and 
Mg II and possibly other metals in WASP-12b bv IPossati et al.l koi(h . 



Of particular importance to o ur study in this paper are the detections of H I made 
for HD 209458b and HD 189733b. Ividal-Madiar et al.l J2003h observed HD 209458b with 
Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST) 
and found a Ly« absorption level of 15 ± 4%. They calculated the Roche lobe of the planet 
to be at 2.7 Jupiter radii, but the height of the absorbing hydrogen required to produce 15% 
absorption to be 4.3 Jupiter radii. This indicate s that the hydrogen overfills the Roche lobe 
and thus is escaping. IVidal-Madjar et al.l (120031 ) also calculated a simple escape model and 



placed a lower limit on the escaping material at 10 g s , the rate is higher if the absorption 
is saturated. 



Additional observations of HD 209458b were presented in IVidal-Madjar et al.l ( 120041 ). 
They found absorption in H I of 5 ± 2%, O I of 13 ± 4.5%, and C II 7.5 ± 3.5%. The 15 ± 4% 



^Note, however, that there is significant dispute over certain Spitzer resu lts, most notably the lTinetti et al 



([20071) detection of H2O in HD 209458b; e.g., see lEhrenreich et al.l (|2007l ) for an alternative analysis of the 
data. 
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H I absorption measurement of IVidal-Madjar et al.l (I20031) is at higher resolution and de- 
tected at significant levels over a limited velocity range; IVidal-Madjar et al.1 ( 20041 ) estimate 
that this translates to 5. 7 ± 1.5% over the entire Lya line at the lower resolution of the latter 
observations, and thus the two measurements are consistent. The O I and C II transitions 
are multiplets including excited states that are ble nded in the lower-resolution data. The 
higher- resolution data of IVidal-Madjar et al.l (120031 ) show that the observed emission from 
the ground state O I line is much weaker than the observed emission from O I* or O I** due 
to some combination of the inher ent stellar emission and int erstellar absorption. This argues 
that the blended emission line in lVidal-Madjar et al.l (l2004f ) is dominated by O I* and O I** 
emission, and the 13% absorption from the planetary atmosphere must primarily come from 
the excited transitions. Similar arguments leave the carbon absorption ambiguous between 
ground-state C II and C II*. Both results imply mini mum densities in th e atmo sphere, be- 
cause these excited states are collisionally populated. IVidal-Madjar et al.l ( l2003l ) argue that 
hydrodynamical escape or "blow-off is the likely mechanism, in which abundant H I escapes 
at a high rate and other species (in this case the O I and C II are dr agged hydrodynamically 
to fill the Roche lobe). However, note that lEkenback et al.l ( l2010l ) offer an alternative in- 
terpretation that involves an inflated hydrogen atmosphere and energetic neutrals (created 
through charge-exchange reactions between the stellar wind and the planet's atmosphere), 
but not evaporation. 



Lecaveher Pes Etangs et al.l (120 lOf ) used the spectroscopic mode of the Advanced Cam- 
era for Surveys (ACS) onboard HST to observe Lya absorption in HD 189733b. They mea- 
sured a Lya absorption level of 5.05 ± 0.75% for in this target; the broadband transit depth 
i s 2.4% , meaning that the differential absorption level is 2.65%. iLecavelier Pes Etangs et al. 
( I2OIOI ) argue that the H I they detect must be unbound, and calculate constraints on the 
escape rate and EUV flux. They flnd a range of 10^ g s~^ to 10^^ g s~^ for the escape rate 
and 10 to 40 times the solar ionizing flux, with the best flt at dM/dt = lO^'' g s~^ and 
F{EUV) = 20Fq{EUV). 

It is important to note that all published cases of atmospheric detections in exoplanets 
are based on a single absorption line (or close doublet) or a single photometric detection (with 
a speciflc species only inferred as responsible) rather than a s uite o f lines, with one partial 
exception: using low-resolution ACS spectra, iBallester et al.l (120071) detecte d absorption in 
HP 209458b corresponding to the Balmer edge of H ll^ lBallester et all J2007h found ~0.03% 
absorption in addition to the planetary disk's baseline absorption. Note, however, that this 



^By the Balmer edge of H I we mean the combination of the n = 2 ^ 00 transition edge at 3646 A, the 
associated bound- free absorption from the n = 2 state at nearby but shorter wavelengths, and bound-bound 
absorption from n = 2— >nS>2at nearby but longer wavelengths. 
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is a blended, conglomerate feature (especially when observed in low-resolution spectra) and 
not a single, resolved line. 

Atomic hydrogen has been detected through its Lya transition in the aforementioned 
targets and in the Balmer edge of HD 209458b, but not in the second-most prominent 
hydrogen line, H«. While the detection of a single line such as Lya is important in its own 
right, it is very limited because it only reveals the presence of that atomic or molecular species 
and an extremely crude measure of density. However, the detection of an additional line 
would provide interesting physical constraints on the density and temperature. For example, 
measuring both Lya and Ha absorption (and assuming an absorption- density relationship to 
derive the density) would enable the calculation of an excitation temperature, Texc- We note 
that Texc does not necessarily equate to the local gas temperature. However, it is still an 
interesting diagnostic of the physical conditions of an exosphere because a difference between 
the local gas temperature and Texc implies a deviation from LTE. 



Winn et al.l ( l2004j ) searched for Ha in HD 209458b with the Subaru High Dispersion 



Spectrograph. They were unable to detect any Ha absorption, but placed an upper limit of 
0.1% over a 5.1 A band (to match the range over which Lya is detected in HD 209458b). In 



turn, 



Winn et al.l (120041 ) placed a limit on the column density of hydrogen in the first excited 
state (n = 2) of N2 < 10^ c m~^ and an upper limit of T^xc < 8000 K. In their subsequent 



study, iBallester et al.l (120071 ) argue that the Balmer edge detection in HD 209458b implies an 
absorption level of 0.3% absorptio n in Ha, though not necessarily over the same waveband 
as measured by lWinn et al.l (120041 ). 



In this paper we report on our search for Ha absorption in f our exoplaneta r y sys- 
tem s previously examin ed for sodium and potassium absorption in iRedfield et al.l (120081 ) 
and I Jensen et al.l (1201 ll ). hereafter R2008 and J2011, respectively. Notably, we detect what 
appears to be Ha absorption in the atmosphere of HD 189733b. This is the first such detec- 
tion of Ha absorption in any exoplanetary atmosphere (or exosphere) and allows us to make 
a unique determination of Texc for the hydrogen associated with this planet. We describe the 
observations and basic data reduction in ^ In [|3]we describe how we analyzed the reduced 
spectra including error determinations. In §l]we describe our basic results. We discuss these 
results in ^and summarize in ^ 
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2. OBSERVATIONS AND DATA REDUCTION 

2.1. Description of Observations 

We obtained observations of the planet-bearing systems HD 147506, HD 149026, HD 
189733, and HD 209458 using tlie Hobby-Eberly Telescope (HET) from August 2006 to 
November 2008. These data were used to make the first ground-based detection of an 
exoplanetary atmosphere by measuring Na I in HD 189733b (R2008). These results were later 
confirmed in an semi- independent analysis of the same dataset (J2011), who also examined 
the spectra for K I absorption. Recent HST observations have also confirmed the presence 
of Na I in HD 189733b (Huitson et al., 2011, in preparation). 

Details of the observations can be found in R2008 and J2011, but here we provide 
an overview. Observations were made with the High- Resolution Spectrograph (HRS) at a 
nominal resolution of i? ~ 60, 000. Individual exposures were limited to 10 minutes in length 
in order to avoid saturation (and are in some cases shorter due to observing constraints). 
Approximately 20% of the exposures for each target were taken while its respective planet 
was in transit. Specific details are provided in Tables [T] (observational statistics) and [2] 
(system and orbital parameters). 

Tabled] also lists the telluric standard stars used for each observation. Telluric standard 
observations were made either immediately before or after primary science observations, thus 
limiting spatial and temporal variations in air mass and water vapor levels. Examples of the 
primary target and telluric standard spectra are shown in Fig. [T] 



2.2. Data Reduction 

The HRS instrument on the HET is an echelle spectrograph. We used standard IRAF 
procedures to remove the bias level and scattered light, flatten the field, extract the apertures, 
and determine the wavelength calibration from the ThAr lamp exposures. We extracted to 
one- dimensional spectra the primary science observations, the telluric standard star obser- 
vations, and the ThAr exposures. 

In order to perform a proper telluric correction, we must remove the imprint of the strong 
Ha line in the telluric standard star. We model a high-order spline (see Fig. EJ through the 
entire spectral order containing the Ha line in all of the telluric observations, then calculate 
an average model to be removed, via division, from these observations. When dividing this 
model from each telluric observation, we also add a low-order polynomial to account for 
wavelength shifts and differing scales between observations. After this division, we are left 
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with a normalized telluric spectrum, free of the telluric standard's Ha imprint. The primary 
science observations are then corrected for telluric absorption using the normalized telluric 
spectra and standard IRAF procedures. We note that we find no apparent telluric residuals 
in the resulting transmission spectra (see §3.ip . 

We separate the exposures by whether they are in- or out-of-transit, and coadd spectra 
from a given group, leaving us with a master in-transit spectrum and a master out-of- 
transit spectrum. Two additional issues arise in the coadding process. First, by measuring 
unsaturated lines in the ThAr lamp exposures, we note that the resolution varies slightly as 
a function of time. Secondly, we must apply small wavelength corrections, both an overall 
shift and a higher-order correction based on cross-correlating the stellar lines. These two 
steps were described in more detail in J2011. The final coadding is weighted by the errors 
at each point in the spectrum, with the errors based on photon statistics and read noise. 



3. ANALYSIS 

3.1. Transmission Spectra 

With our master in- and out-of-transit spectra we derive a transmission spectrum. We 
perform a wavelength correction similar to that used in our coadding process (an overall 
shift plus a higher-order correction from cross-correlating stellar lines). The spectrum is 
normalized, with the continuum defined several angstroms away from the line center. We 
then calculate the transmission spectrum. 



We calculate a measurement of the absorption, Mabs, which we define as 

where the c, b, and r subscripts represent the central primary integration region, and the 
blue and red comparison integration regions, respectively. 



3.2. Estimating Errors with an Empirical Monte-Carlo Method 

Because of the many steps in our reduction process, we suspect that systematic errors 
may be larger than the formally propagated statistical errors. These systematic errors may 
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come from either our reduction process or from astrophysical sources. Systematic errors from 
our reduction process would include the telluric correction step, normalization of spectra, and 
our wavelength corrections. The dominant astrophysical source of systematic error would 
be stellar variability. To quantify the effect of these systematic errors, we use an "Empirical 
Monte-Carlo" method (hereafter EMC). The details of this method can be found in R2008 
and J2011, but we summarize here. We use our coadding process described in §2.2l on various 
subsets of the in- and out-of-transit data (always of the same type) and compare them to the 
master in- or out-of-transit spectra using our procedure for calculating transmission spectra. 

The three important combinations we examine are: (1) in-transit subsets vs. the out-of- 
transit master spectrum or "in-out" ; (2) in-transit subsets vs. the in-transit master spectrum 
or "in-in"; and (3) out-of-transit subsets vs. the out-of-transit master spectrum or "out-out." 
Each method will have some characteristic distribution, but we expect that the "in-out" 
method will be centered around the master transmission spectrum measurement and the 
"in-in" and "out-out" methods will be centered around zero. Additionally, we expect that 
the distribution of the "out-out" method will provide the best estimate of the true systematic 
error (see R2008 and J2011 for more details). 



3.3. Confirming Detections with a Control Line 

Another important check on our results is to see whether they are the systematic result 
of our coadding and transmission spectrum calculations near a strong stellar line. To check 
this, we performed the above procedures on a strong stellar line of Ca I at 6122 A, which 
we do not expect to manifest in exoplanetary atmospheres (calcium should condense out 
at the temperatures implied by models, e.g., into CaTiOs). If our method is robust even 
when performed on such a line, the results of the master measurement and the EMC distri- 
bution should be consistent with zero (see R2008 and J2011 for additional details). These 
calculations were performed by J2011, and we do not repeat them here (see Fig. 2 and 5 of 
J2011). In J2011, we did indeed find the master measurements and EMC distributions to be 
consistent with zero for all four systems, an indication that any absorption results on other 
lines are not simply an inevitable, or even likely, consequence of our analysis procedures. 



4. RESULTS 

The resulting transmission spectra from the reduction and analysis processes described 
above are shown in Fig. [31 The corresponding EMC analyses are shown in Fig. HJ and include 



■'in-out," and "out-out" results. In both figures, each panel corresponds to one of 
our targets. For the commensurate results of the control line analysis of Ca I described in 
§3.3^ see Fig. 2 and 5, respectively, of J2011. 

Initially, we performed our integration (i.e., the calculation oi Mabs from Eq. |2]) over a 16 
A band. We anticipate that Ha may be broad, and therefore chose this in tegration band to be 



broad er than the 12 A band used to measure Na I in R2008, J2011, and ICharbonneau et al. 



( I2OO2I ) and the 8 A bands used to measure Ca I and K I in J2011. However, because we 
compare this measurement to adjacent bands of equal width, a 16 A integration band requires 
a continuous spectrum of 48 A centered on the line. The spectral width of the echelle orders 
is ~100 A, and our selection of integration bandwidths are limited by this width and how 
close the line is to the center of the order. However, we also define different bandpasses 
and perform the corresponding measurements for both the master measurement and the 
measurements of the EMC calculations. In Table [3] we show all the resulting Mats for a 16 A 
bandpass, which is the bandpass that is represented in the master transmission spectra and 
EMC analyses that are shown in Fig. |3]-|11 For comparison. Table |3] also shows the results of 
a 6 A bandpass, which is discussed in §4.1[ A few additional custom bandpasses are defined 
to examine the unique feature in HD 209458b, which we describe in §4.31 

We also discuss stellar variability as it pertains to the features noted in HD 189733b 
and HD 209458b. Our observations cover a wide range in time and allow for period analysis 
at the level of individual observations. We find significant stellar variability at H« in the 
case of the star HD 189733, and this variability must be understood in order to properly 
interpret the transmission spectrum of the planet. 



4.1. HD 147506b and HD 149026b 

The transmission spectra of both HD 147506b and HD 149026b are relatively featureless 
compared to the noise level. Notably, the strongest "noise" feature is at approximately 6553 
A in HD 147506b, which, at our settings, corresponds to the known location of a large 
group of bad pixels (roughly a 3 x 3 square of pixels) on the CCD chip. Standard IRAF 
procedures and the telluric correction step largely take out the effects of these bad pixels in 
most individual exposures, but inspection here of the final transmission spectra shows that 
some residuals remain that are only obvious at high precision. 

Though the transmission spectra for these two targets do not show any clear features 
that appear to be absorption or emission, the master integration measurements and corre- 
sponding EMC results show a slightly more complicated situation. The master integration 
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ineasureinent for HD 147506b results in an absorption level that is 2. 80" significant relative 
to the width of the "out-out" EMC, which we assume is the best estimation of the total 
systematic and statistical error (see §3.21 and similar discussions in R2008 and J2011). No- 
tably, this master measurement is not well-aligned with the "in-out" distribution, which is 
an exception in the 16 transmission spectra calculated between this work and J2011. The 
best-fit centroid of the "in-out" distribution is only 0.7a away from zero, again using the 
Gaussian a from the "out-out" distribution as the error. The "in-in" distribution is also not 
quite centered at zero (though it is only approximately l.Ocr away from zero using its own 
Gaussian width). 

There are a few things we can say about the results for HD 147506b. First, based on 
the master transmission spectra, it appears to be the only one of the four targets where the 
aforementioned section of bad pixels, affecting data points near ~6553 A, comes into play 
in the final transmission spectrum. Second, as we noted in J2011, HD 147506 is the faintest 
star in our sample, with the resulting lowest S/N (comparable to observations of HD 149026 
but much lower than HD 189733 and HD 209458). It also has the broadest and shallowest 
stellar features of the four s tars in our sample, due to its v sin i of 20.8 ± 0.3 km s~^, which is 



> 3x the other three stars (jPal et al.ll2010l : lButler et al.ll2006l ). This complicates our process 



of cross-correlating stellar lines of different spectra in order to make additional corrections 
to the wavelength solution (note, for example, that there is also a small alignment artifact at 
~6575 A in this target). Therefore, it would have been reasonable to identify this target a 
priori as the most qualitatively uncertain, a conclusion we also reached in J2011 where similar 
issues were seen in the transmission spectrum of this target. Finally, we note that we obtain 
different result when we integrate over a narrower bandpass. Our default bandpass is 16 A; 
we also calculate a 6 A bandpass, which is chosen so that the blue comparison band misses 
the ~6553 A noise feature. This results in the master transmission spectrum integration 
being reduced; note that this is an average value, so the two are already normalized despite 
the different bandpasses. The distribution of the "out-out" EMC is also wider, so the 
significance drops to 0.58cr. The consistency of the master measurement and the "in-out" 
distribution is slightly improved in an absolute sense, and dramatically improved relative to 
the now-broader widths of the both the "in-out" and "out-out" distributions. Furthermore, 
note that in both cases, the ~6553 A noise feature may affect the overall normalization of 
the transmission spectrum. 

In the transmission spectrum HD 149026b there is a slight "emission" bump just to 
the red of the stellar Ha line center. This turns out to be significant at the 2. Oct level. 
Unlike the unusual EMC distributions (in that there are apparent internal inconsistencies) 
of HD 147506b described above, in this case the "in-in" and "out-out" distributions are both 
consistent with zero, and the "in-out" distribution is consistent with the master integration 



-10- 



ineasureinent. Changing the bandpass also does not affect these results in the same way that 
we described above for HD 147506b. The feature is visually subtle and only 2. Oct, so we are 
disinclined to associate this with a real feature. However, if the feature is real, the primary 
way that an emission feature could be observed in a transmission spectrum is by induced 
stellar activity from a star-planet interaction (SPI), which might occur on the portion of the 
stellar surface facing the planet. We discuss this more in §5. 2. H and §5.3.31 



4.2. HD 189733b 

We detect a dramatic absorption-like feature in the transmission spectrum of HD 189733b, 
with a value of Mabs = (1-06 ± 0.24) x 10~^ (integrated over 16 A and Mats defined as in 
Eq. |2] and using the standard deviation of the "out-out" EMC distribution as our error). 
This detection is significant at a 4.4o" level. The EMC results place important constraints on 
any otherwise undetected systematic errors, including any remaining atmospheric residuals 
of our telluric correction and planet-independent stellar variability. The statistical signifi- 
cance of this result argues strongly that this is a real, transit-dependent feature rather than 
a reduction or analysis artifact. We also note that when we inspect the transmission spec- 
tra of the several thousand EMC iterations visually, the results are qualitatively consistent 
with the measurements: all "in-out" iterations consistently show an absorption feature in 
transmission, whereas the "in-in" and "out-out" iterations show variation in the residuals of 
emission or absorption in the Ha line, but the variation is centered about zero and smaller 
in magnitude than average magnitude of the "in-out" absorption. The variation in the "in- 
in" and "out-out" iterations is probably due to stellar variability, which we discuss below. 
Another qualitative argument in favor of the HD 189733b result being a real phenomenon 
associated with transit is the null results in HD 147506b and HD 149026b. Similar to the 
Ca I control line in these targets, this shows that an artificial feature is not an inevitable 
result due to some effect of our transmission spectrum procedure performed on a strong 
stellar line. 

However, the level of variability in the "in-in" and "out-out" — which we noted above 
was less than the magnitude of the absorption feature in the master transmission spectrum — 
is worth further inspection. The EMC takes a random sample of the observations and is 
unlikely to sample all observations from a given night. Therefore, if there is significant stellar 
variability on time scales greater than the length of the individual observations but shorter 
than the time spanned by different sets of observations (e.g., nights), then the EMC as 
initially constructed may not adequately characterize the variability. Therefore, we make an 
additional attempt to do so by calculating a transmission spectrum and integrated absorption 



- 11 - 



measurement for each individual out-of-transit observation relative to the master out-of- 
transit spectrum. This assumes that the master out-of-transit spectrum sufficiently samples 
any variability and averages it out; using only out-of-transit observations explicitly removes 
the dependence of planetary transit from the calculation. This provides us with a time series 
representation of the variability, on which we can perform period analysis. 

We used a Lomb-Scargle periodogram to analyze this time series, which is shown in 
Fig. O along with the window function associated with our data. There are many periods 
with power spectral density (PSD) of greater than 3a significance, according to both analytic 
determination and white noise simulation estimations of the false alarm probability. The 
strongest group of peaks in the periodogram, in terms of PSD, is at high frequency (~1 day); 
the second strongest set of peaks is near ~12 days; the third strongest set of peaks is at ~8.5 
days; and more are found at one-half of the ~12 day group, as well as lower frequencies of 
>16 days. In contrast, similar periodograms for HD 147506 and HD 149026 do not show any 
variability that exceeds the 3a false alarm threshold (we deal with HD 209458 in §4.3p . Given 
that HD 189733 is the most chromospherically active star in our small sample (see discussion 
in §5.2.5p . this might be expected. Of th e periods just described, the ~12 day period is the 
one with a clear physical interpretation. iBouchy et al.l ( l2005l ) estimated a rotational period 
of ~11 days, based on the v sini and the assumption that we are vie wing the star along its 
plane of rotation. Similarly, iHebrard fc Lecavelier Des Etangsl (120061 ) found microvariability 
at several periods, with the dominant period to be at 11.8 days, whi ch they interpreted as the 
rotation period. "Quas i-periodic" flux va riations were measured by lWinn et al.l (l2007al ) with 
a period of 13.4 days. iFares et al.l (120 lOl ) suggest the discrepancy between these latter two 
results may be due to the temporal variation in stellar spots coupled to differential rotation, 
and flnd an equatorial period of 11.94 ± 0.16 days. 



Fares et al.l (120 lOl ) observed Ca II and Ha and found signiflcant variability. For each 



line, they observe over three epochs and calculate a best-flt period. The Ha best-flt periods 
are 11.6+01, 14.4+i:^, and 13.8?:|3 days for observations in 2006 June and August, 2007 
June, and 2008 July (though the authors note that the phase coverage is poor for the 2006 
data). Note that all three best-flt Ha periods are consistent with our period of 12.09 days 
(the period with the largest PSD between 10 and 15 days) to within 2a; 12.077 days is also 
consistent to within la with their derived equatorial period. 

Notably, there is also a straightforward explanation for the high-frequency peaks. The 
window function peaks at 0.9967 days, which is expected because the timing of the HET 
observations is correlated with the length of the sidereal day, which is 0.9972 civil days. The 
two peaks in the periodogram that are closest to this peak in the window function are also 
the two most powerful peaks in the periodogram, at 1.086 and 0.920 days. These peaks are. 
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to high precision, plus and minus one 12.077 day period in frequency space away from the 
window function peak. This strongly suggests that these peaks are due to aliasing of the 
longer period in conjunction with the window function peak of 1 sidereal day, and therefore 
the power in these peaks should be disregarded. 

We attempt to remove the variability from our master transmission measurement by 
making a fit based on only the out-of-transit variability. Then we can estimate an additive 
calibration term for any observation (in- or out-of-transit) by phasing it to this function. 
These calibration terms are then averaged for all in-transit observations in order to calculate 
a calibrated master transmission spectrum integration, and also for each iteration of the 
EMC. We show a fit to data phased to the 12.077 day period in Fig. [6l using a low-order 
Fourier series. In Fig. [7] we show this same fit as a function of absolute time. Both Fig. [H] 
and [7] show the general tendency of the in-transit residuals to be below the out-of-transit 
residuals. When we calculate a new periodogram with the residuals, the PSD of the ~ 1 
day high-frequency peaks is greatly reduced — the two highest-power peaks disappear and 
the highest remaining ~ 1 day period with the largest PSD is at the peak of the window 
function. The other significant peaks of higher frequency than the 12.077 day period that we 
noted above are also significantly reduced. We limit the fit to two terms in the Fourier series 
(i.e., the function is a linear combination of sines and cosines with periods of P and P/2 
where P is the 12.077 day period) based on what we see in the periodogram; if more terms 
are allowed, the function begins to fit for some of the periodicity caused by our sampling 
frequency as characterized by the window function. 

The calibrated EMC and master measurement gives us a new result of (8.72 ± 1.48) x 
10"^. The error from the EMC in this case has gone down significantly because we are 
explicitly fitting for the dominant source of astrophysical noise. HD 189733 is the brightest 
star in our sample with the greatest number of observations and therefore should have 
the highest S/N and lowest corresponding error if astrophysical noise is removed. Prior to 
calibrating the error, it was largest for HD 189733b. The calibrated error is smaller than 
for HD 147506b and HD 149026b with their fainter central stars, but larger than for HD 
209458b and its comparably bright star. Using the uncalibrated measurements, our result 
was significant at a 4.4o" level; the calibrated result is significant to 5.9a, and the EMC 
results behave as we expect — the widths of the three EMC distributions are each reduced by 
a reasonably consistent fraction of between 39 and 46%, the master "in-out" measurement 
is consistent with the "in-out" EMC, and the "in-in" and '"out-out" EMCs are consistent 
with zero. 
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4.3. HD 209458b 

HD 209458 is equally intriguing for different reasons. We do not detect an overall 
absorption or "emission" signal, as the master integration measurement over the 16 A band 
is consistent with zero. All three variations of the EMCs are centered on zero, although 
the "in-in" and "in-out" distributions are slightly asymmetric. However, the transmission 
spectrum shows a dramatic feature with a "spike" to the blue of stellar line center and a 
"dip" to the red. Both features peak at roughly 0.5% absolute deviation from zero, and are 
several angstroms wide. The feature is roughly symmetric if a reflection is applied about the 
zero point and line center. When we use our EMC process and shift the integration bands to 
correspond to these features, the EMCs confirm that these are consistent features that are 
highly significant. In order to explore these features which appear to reflect about the stellar 
Ha line center, we examined <St> in two 5 A bands with their edges centered on the stellar 
line center (unlike the 16 A and 6 A bandpasses of Table [3] which are centered on the stellar 
line); 5 A is the approximate width of the spike and dip. We find that in the blueward 
spike, <St> (omitting the comparison bands) is an 6.6cr "emission" feature, and in the 
redward dip, <St> is a 6. la absorption feature. The corresponding "in-in" and "out-out" 
distributions are consistent with zero to within better than 0.5o". Visual inspection of the 
EMC iterations qualitatively confirms these results — i.e., these features are consistently seen 
in the various "in-out" iterations and not present in the "in-in" and "out-out" iterations. We 
also note that if we i ntegra te over the range of —130 to -1-100 km s~^ — the range over which 



Vidal-Madiar et al.l (120031 ) detected significant Lya absorption, and the range over which 
Winn et al.l (J2004f ) search for Ha absorption — we arrive at a value for <St> of 1.27 x 10^^. 



The positive value is due to the asymmetry of the velocity limits favoring the blueward 



spike; over a similar range, I Winn et al.l fl2004f ) placed an upper limit of 0.1%. We suggest an 



explanation for the discrepancy later in this section. 

It is difficult to explain apparent "emission" features in a transmission spectrum. How- 
ever, it is not impossible for such features to occur. Because of the nature of normalizing 
St to 0, an "emission" feature indicates a wavelength-dependent apparent planetary radius 
{R\/Rp) of less than one. However, this assumes a constant stellar disk intensity. If the stel- 
lar disk intensity changes as a function of both wavelength and time, wavelength-dependent 
features may be observed in the cumulative transmission spectrum. For example, in the 
Rossiter-McLaughlin effect, the exoplanetary disk blocks different portions of the stellar 
disk, and each portion of block stellar surface has a different radial velocity relative to the 
observer due to stellar rotation. This creates wavelength-dependent variations in the disk 
intensity. The apparent symmetry of the HD 209458b feature also indicates the possibility 
that it is due to some sort of velocity misalignment between two absorption features that 
is not fully corrected in our transmission spectrum calculation (cf., the artifact at ~6575 A 
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in the transmission spectrum of HD 147506b). In the case of HD 209458b, this would have 
be a very broad feature that is misahgned between the in- and out-of-transit spectra. The 
misahgnment would also need to be relative to the core of the stellar Ha line, as there is no 
residual indicating that the strong central Ha cores are misaligned in our spectra. While we 
cannot completely rule out a reduction artifact, if such a misalignment is responsible for this 
feature, its consistency in the various "in-out" EMC iterations argues for a genuine physical 
origin. This could perhaps be due to SPI-related stellar variability. Another hypothesis is 
that the "broad" component of the stellar Ha line actually includes absorption from excited 
hydrogen that has been lost by the planet, remains around the star in some sort of toroidal 
shape, and is influenced by the planet's orbit. We discuss potential interpretations more in 
^1531 



As with HD 189733b, we investigated whether or not there is variability beyond just 
transit that is not picked up by the EMCs. This is represented in Fig. [8], using the 5 A 
integration band to the blue of the Ha line as the measure of Mabs- Positive measurements 
therefore mean that there is a blue spike and a red dip relative to Ha; negative measurements 
would indicate an inverted feature. Notably, there does seem to be a phase correlation, as 
the blue spike is prominent during transit and near secondary eclipse (orbital phase ~ 0.5), 
but at intermediate values (0 ~ 0.25 and (j) ~ 0.75), there is either no feature or an inverted 
feature (we examined the spectra and confirmed visually that these descriptions are the 
correct interpretation of Mabs, which again, in this case, refers to the offset measurement 
centered on the blue spike). Such features apparently do not show up as readily in the "in- in" 
and "out-out" EMCs due to the fact that observations are selected individually and there 
are very few, if any, iterations that are dominated by a single night. If the phase correlation 
seen in Fig. [8] is real, and the effect seen in the master transmission spectrum is not strictly 



a tran sit effect but a broader phase effect, then it also explains the fact that IWinn et al. 



(120041 ) did not note this feature, as they made observations over only ~ 5 hours on one night, 
through a single transit, and did not have the wide range of phase coverage that we present 
here. 



5. DISCUSSION 
5.1. Limits on Ha in HD 147506b and HD 149026b 



As described in §4.H in neither HD 147506b nor HD 149026b do we detect clear Ha 
absorption. The atmospheres of these planets are not well-studied, with no Lya detections 
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(or null results) jj and comparatively little IR work. No Spitzer observations have been 
taken of HD 147506b f Seager fc Demingj 12010), and only 8 /im Spitzer observations have 



been taken of HD 149026b (JHarrington et al.l 120071 : iKnutson et al. 



20091 ). In an analysis of 



the potential connectio n between stellar act ivity and the presence or lack of atmospheric 
temperature inversions, iKnutson et al.l ( l2010l ) categorized these two planets' atmospheres as 
"unknown" with respect to temperature inversions. 

In J2011 we examined these spectra for absorption at Na I AA5889,5895 and K I A7698. 
We did not find clear evidence of absorption from either species, though we reported the 
suggestion of Na I absorption in HD 149026b, significant at the 1.7cr level, with a self- 
consistent EMC and plausible features in the transmission spectrum. The measured strength 
of this possible absorption would very tentatively predict no temperature inversion, as the 
extra opacity source(s) needed to produce temperature inversions may tend to suppress the 
observabihty of atmospheric features, e . g., th e suppressed Na I absorption in HD 209458b 
Jcharbonneau et al.l[2002l : ISnellen et all2008h compared to HD 189733b (R2008; Huitson et 
al., 2011, in preparation). 

The lack of Ha detections in these two targets as well as lack of comparison points in their 
atmospheric spectra in the literature of these planets severely limits our ability to draw firm 
conclusions about their atmospheres. In §5.2[ we discuss limits on the excitation temperature, 
Tf^xci of hydrogen due to the now-available combination of Lya and Ha measurements for HD 
189733b and HD 209458b. With Lya detections, we could in principle put similar constraints 
on Tf,xc in these two targets. 

We also note here that the a priori expectation of detectable atmospheric absorption 
in HD 147506b is by far the lowest, and also lower for HD 149026b than the remaining 
two targets. This is primarily based on S/N limitations due to the fact that both central 
stars are fainter than HD 189733 and HD 209458, which have similar apparent magnitudes — 
HD 149026 is fainter by ~0.5 mag and HD 147506 by ~1 mag. In addition, HD 147506b 
also has by far the smallest scale height of our four targets due to its very high mass even 
though its radius is comparable to the other planets. In terms of scale height, HD 149026b is 
comparable to HD 189733b, but the prospect of detecting absorption is made more difficult 
due to the lower S/N resulting from the respective apparent brightness of their central stars. 
Despite these relative limitations, these systems have bright enough central stars that if 
Lya absorption exists in their exospheres, it should be detectable with a reasonable amount 



■^Thcre are observations of HD 149026b using ACS in the HST archives that cover Lya. However, these 
observations were taken while STIS — which is, generahy speaking, preferable to ACS for spectroscopy and 
therefore for detecting wavelength-dependent absorption in the UV — was not operating, and there are no 
published results associated with this program. 
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of observing time. Furthermore, we note that our upper hmits on Ha are a factor of a few 
smaller with respect to the HD 189733b detection and the HD 209458b feature. HD 147506b's 
scale height is much smaller but its Roche lobe is much larger than the other targets; this 
reduces the possibility that a true, unbound exosphere exists. Our non-detection of Ha 
does not confirm this, but is certainly consistent with it, and therefore not unexpected. HD 
149026b is different in that its scale height and Roche lobe are comparable to HD 189733b 
(slightly larger in both cases). Therefore, our non-detection of Ha in this target is arguably 
indicative of a more interesting difference between this target and HD 189733b and HD 
209458b. Whatever processes that create the features seen in these two targets (see §5.21 and 
§5.3p are apparently not present, or significantly weaker, in HD 149026b. 



5.2. The HD 189733b Ha Detection 

5.2.1. Identification and Description 

We argued in §4.21 that the feature we observe in HD 189733b is highly correlated with 
planetary transit and is a real feature of the transmission spectrum. If this premise is 
correct, two basic possibilities exist: that the feature is the result of the influence on the 
stellar Ha line from SPI, or that it is the result of absorption from the planet's atmosphere or 
exosphere. HP 1897 33b is the most chromospherically active of the stars in our sample (see 



Knutson et al.ll2010l . and also the discussion in §5.2.51) . and thus the first possibility should 
not be ignored — especially given that we strongly consider SPI as a possible explanation 
for the unusual feature in HD 209458b in §5.3[ However, confirmed instances of SPI are 



quite rare (e.g., iPoppenhaeger et al.ll201ll . and references therein). Furthermore, in cases 
where SPI has been detected, the effect on the star is complex as a function of time. For 
example, HD 179949, a non-transiting system, shows the influence of the planet for nearly 
half of its orbital cycle, a nd also may show ad ditional variability on larger time scales such as 



the star's activity cycle fjShkolnik et al.ll2008l ). Recently, the possibly of HD 189733' s X-ra y 



activity being correlated with the planetary phase was suggested by iPillitteri et al.l (120 111 ) 



who acknowledge that further observations are needed to confirm or reject this hypothesis. 

In the simplest case, though, we might expect that the planet influences additional 
activity on the portion of the stellar surface that is facing the planet. If this were the case in 
HD 189733, the in-transit observation would have additional chromospheric emission during 
transit that "fills in" the Ha absorption core. This would result in an emission feature 
in the transmission spectrum rather than the absorption feature we see. Therefore, the 
interpretation that the absorption comes from some portion of the planet's atmosphere is 
highly preferred. 
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As discussed in ^ there are many known detections of exoplanetary atmospheres, and 
this includes exospheric absorption from the Lya transition of hydrogen in HD 189733b and 
HD 209458b and the Balmer edge of H I in HD 209458b. The detection of Ha in HD 189733b 
is the first such detection of this hne in an exoplanetary atmosphere or exosphere. It is also 
more generally the first detection of atomic material in an exo planetary atmosphere that is 



i n an excited state via a resolved line, and combined with the iLecavelier Des Etangs et al. 



( I2OIOI ) detection of Lya, results in the first case where a set of multiple, resolved lines (which 
do not come from the same multiplet or an ambiguous absorption edge) have been observed 
in a single exoplanetary atmosphere. 

We detect significant Ha absorption in HD 189733b over a range of —24.1 to 26.6 km s^ . 
This is determined by the continuous range over which St < —3a, with the a here being the 
standard deviation of the nearby spectrum. We r iote that this is mu c h nar rower than the 
velocity range of —13 to 100 km s~ over which IVidal-Madjar et al.l ( 20031 ) detected Lya 
in HD 209458b. The Ihecavelier Des Etangs et all fcoioh detection of Lya in HD 189733b 
used ACS at a resolution of i? ~ 100. Therefore, the Lya line was not resolved, and we are 
unable to ascertain how our detection compares to Lya in this target (though we discuss 
this more below). However, the contrast with the IVidal-Madjar et al.l (120031 ) detection is 
still intriguing, because it implies at least one of the following: (1) the Lya absorption of 
these two targets have very different velocity ranges, or that (2) Lya and Ha absorption in 
HD 189733b are not over the same velocity range. 



5.2.2. Velocity Range 

The narrow velocity range over which we detect this absorption is also interesting be- 
cause, at least superficially, it raises the question of whether or not we are detecting hydrogen 
in the exosphere. With R^d i89733b = 1.14/2jupiter and Mhd i89733b = 1.138Mjupiter, the escape 
velocity. 



2GM 
R 



(3) 



is 59.5 km s~^ (at the planetary surface). With the ILecavelier Des Etangs et al.l (120 lOh detec- 
tion of Lya, the absorption was not resolved over the 150 km s~^ full-width at half - maxim um 
(FWHM) of the Lya line in the ACS spectrumo ILecavelier Des Etangs et al.l ( 2010[ ) ar- 
gue that under any scenario, the hydrogen must be escaping — though unresolved, the total 



"^The instrumental resolution is ^300 km s ^ (FWH M) , but Lecavelier Des Etangs et al.l ( 20101) estimate 
the stellar Lya width from the empirical relationship of [Landsman fc SimonI (|1993I ). 
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aniount of absorption demands that either some of the gas is at a velocity greater than the 
escape velocity, or some of it is beyond the planet's Roche lobe. 

Our detection of Ha is not as clear in this respect. The maximum absolute velocity 
at which we detect clearly Ha absorption, ~25 km s~ , is significantly smaller than the 
escape velocity at the planetary surface of ~59.5 km s~^. The escape velocity is equivalent 
to 25 km s~^ at approximately 5.7-Rhd i89733b, or 6.4i?jupitor (which would also be outside 
HD 189733b's Roche lobe). The maximum apparent planetary radius implied by the Ha 
absorption is 1.35i?HD i89733tlj; which is not outside the Roche lobe. Note, though, that 
calculating this apparent planetary radius assumes optically thick absorption; if the Ha 
absorption is optically thin, it could be physically present at much greater radii. Assuming 
the hydrogen responsible for the Ha absorption has the same physical extent as the hydrogen 
responsible for the Lya absorption, our result argues for the former case of lower-velocity 
hydrogen that overfills the Roche lobe rather than gas within the Roche lobe that exists at 
speeds greater than the escape velocity. 



Due to the low resolution of the iLecaveher Pes Etangs et al.l (120101 ) ACS spectra, we 



cannot directly compare the velocity ranges of the Lya and Ha absorption in HD 189733b. 
However, we note that the interstellar Lya absorption is probably satura ted over the velocity 



range of our Ha detection (see Fig. 10 and the corresponding discussion in lLecavelier Pes Etangs et al 



20101 ). Therefore, in order for there to be any detectable Lya absorption in HD 189733b, 
it must indeed cover a wider velocity range than our Ha detection, and any attempt to 
understand how hydrogen in the n = 2 state is formed and sustained in the HD 189733b 
exosphere must recognize that Lya is absorbed by hydrogen over a larger velocity range. 
This also complicates the calculation of column density and excitation temperature that we 
perform below in §5.2.31 and §5.2.41 



5.2.3. Column Density 

If we assume that the Ha absorption is optically thin, we can calculate an estimated 
column density from our absorption measurements by the equation 



\ A J mc^ 



^A I — ) = zrT2h.Xl^N2 , (4) 



^This assumes no calibration as we calculated at the end of tj4.2l Assuming this calibration value is the 
same for the maximum absorption point as it is for the integrated measurement, the maximum apparent 
planetary radius is 1.29i?HD 189733b- 
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where [AA/A] is the fractional area of the stellar disk covered by the absorbing material, 
/23 is the oscillator strength of the Ha transition, and N2 is the column density of hydrogen 
in the n = 2 state. Wx is the equivalent width of the absorption, defined as 



W, 



out -* in 



out 



d\ . (5) 



However, we earlier defined the "transmission spectrum," St, in Eq. [T], such that we can 
substitute and get 

Wx= f-STdX, (6) 



which for approximately equally spaced points in St, and substituting Eq. [21 means that 

Wx ^ -MahsAX , (7) 

Note that the substitution from Eq. [2] is not exact in that Mats only approximates <St>, 
because the former also includes normalization from the comparison bands. In Table HI we 
show the values of Mabs for HD 189733b in three additional integration bands beyond those 
shown in Table [3l Table H] also shows the resulting Wx of these three values of Mabs- As 
should be the case, the Wx is constant to within the errors, because the additional integration 
region of the wider bands is only adding, in principle, <S't>~ 0. 

To calculate N2 is fairly straightforward with the exception of calculating (AA/A), 
because, as we have previously noted, the Lya absorption — which would be our best indi- 
cator of the extent of the hydrogen envelope — is not resolved. Using an equivalent width of 
0.0128 A from the top row of Tab le HI implies a column density of 6.4 x 10^° cm^ (AA/A). 



Lecavelier Pes Etangs et al.l (120101 ) quote values for the absorption, which is equivalent to 
[AA/A), as small as 5% over the entire unresolved Lya line, and as large as 30% if restricted 
to the velocity range of ±49 km s~^ (the escape velocity at 1.65-Rjupitor, the minimum size of 
the hydrogen cloud allowable by the absorption). Thus our effective n = 2 hydrogen column 
density, A''2, could be anywhere from 3.2 x 10^ cm^ to 1.9 x 10^° cm^. Notably, there are still 
strong assumptions in these values. The 30% figure could be even larger if the velocity range 
of Lya is even more restricted than ±49 km s~^, meaning these calculations are upper limits. 
We have also argued that Ha absorption occurs over a smaller velocity range than Lya, and 
therefore may not share its physical extent. This means that {AA/A)-^^ < {AA/A)i^ya, and 
the prior calculations are instead lower limits for a given (AA/A) Lya. 
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as 



5.2.4- Excitation Temperature 
The excitation, or population, temperature of an atomic or molecular species is defined 

AE21 



92 



■exp 



knT, 



B-'- exc 



Solving for Tg; 



Ni gi 
and plugging in constants we obtain the relationship 

l.lSxlO^K 

J. f>, 



In 



4iVi 



(9) 



We have already seen above how many assumptions are in the calculation of A'^2; but we 
can simplify the calculation i n our approximation for A^i . The resolution of the ACS in 
spectroscopic mode as used by iLecavelier Pes Etangs et al.l ( I2OIOI ) is i? ~ 100 (FWHM), or 
~12 A. They detect 5.05% absorption over this line, but if normalized to the baseline transit 
depth of 2.4%, the differential absorption is only 2.65%. Therefore the equivalent width 
is 0.32 A. If the 6-valuqj of the exospheric hydrogen is large enough that the absorption 
is not saturated then we can approximate N1/N2 with M/a, Lya/W^A, Ha- This results in a 
temperature of 2.6 x 10^ K. Of course, we are again violating our previous conclusion that 
Lya and Ha are not fully spatially coincident. We can explicitly state this assumption and 
write the excitation temperature as 

1.18 X 10^ K jj^j 



T 

-J- p' 



4.61 + In 



AA 



Lyg 



AAt 



In this case, even a spatial discrepancy of a factor of 3 between Lya and Ha changes the 
excitation temperature by <25%. 

It is very likely that the Lya absorption is indeed saturated, which would make A^i 
larger and therefore our calculation of T^xc is an upper limit. To quantify this would require 
another additive term of In (A''i^ actuai/A'^i, unsaturated) in the denominator of Eq. [TOl Fig. [9] 
shows the dependence of Tf^^c on A''i (and implicitly on 6- value) as well as A^Lya/AylHa. 
However, there are two important fixed assumptions in Fig. [HI one is that the Lya W\ is 
not corrected for ISM absorption, and the other is that the Ha absorption is unsaturated. 
Both assumptions mean that the N\ and A^2 column densities are lower limits, which would 
respectively decrease or increase the calculated Tf,^c- Notably, iBallester et al.l ( 120071 ) argue 
that any Ha absorption in HD 209458b would be optically thick and in a thin layer (or 
A^Lya/AylHa ^ 1 in our discussion here); if this is true of HD 189733b, these two facts 
would have competing effects on Tga-c relative to our assumptions. 



^Thc term "6-value" refers to the Dopplcr parameter, h = crv2, of a Gaussian velocity distribution. 
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5.2.5. Physical Interpretation 

What sort of physical model can explain these conditions? If we assume that HD 
189733b has zero albedo and the surface is in radiative equilibrium with the star as its only 
energy source, we can calculate a temperature of 1700 K (which we will refer to as Tgg). This 
is an order of magnitude smaller than the excitation temperature that we have calculated. 
The uncertain quantities that would reduce our final calculation oiTexc are A'^i, (AA/y4)Ha, 
and {/S.A/A)i^ya.] T^xc depends on the logarithm of these numbers. Therefore, we would 
need to be off by a total of ~28 orders of magnitude in these quantities in order to make 
Texc = Teg = 1700 K. If we have underestimated the effect of these quantities by a total of 
only five orders of magnitude, T^xc is still greater than T^g by a factor of more than four. 

The discrepancy between Tgg and T^xc argues first that the hydrogen that produces our 
observed Ha absorption must be in an extended, low density portion of the atmosphere. If 
it were too close to the surface, it would likely be in thermal equilibrium with the planet at a 
value close to Tg^; if it were at higher density, it would also likely be collisionally de-excited. 
Therefore, the likely source of excited n = 2 hydrogen is through radiation. Of our four 
systems, HD 189733 is the latest spectral type, an early K star (HD 149026 is a GO star, 
while HD 147506 and HD 209458 are F8 stars). However, it also has the smallest planet- 
star separation, and has the most active Ca II H & K emission. HD 189733's measured 
value of S'hk — a quantity indicating stellar activity — is 0.508, com pared to 0.191, 0.152 , 



and 0.160 in HD 147506, HD 149026, and HD 209458, respectively flKnutson et al.l 120101 ). 
Similarly, HD 189733's value of log (-Rhk) — ^ quantity derived from S'hk that allows for better 
compari son between dif f erent spectral types — is —4.501, compared to —4.780, —5.030, and 



—4.970. iKnutson et al.l (120 lOf ) also note that HD 147506, with the second highest value of 
log (-Rhk) i^ ^^^ target list, is at the edge of the B — V range over which the calibration 
between S'hk and log (-Rhk) is valid. They subsequently argue, based on its spectra at Ca II 
H & K, that HD 147506 has a quieter chromosphere than its derived log (-Rhk) value would 
otherwise suggest. 

The importance of Ca II H & K emission is that it is strongly correlated with general 
chromospheric activity, and we would expect it to therefore correlate strongly with Lya and 
UV activity, although the directly measured evidence for this is limited by a lack of surveys 
in the literature. Stellar Lya and UV activity would be able to populate the n = 2 state of 
hydrogen. Lya emission would directly create n = 2 hydrogen in the planet's atmosphere, 
while higher-frequency UV emission would populate the n > 2 states, which would then 
radiatively cascade down to n = 2. 

If we assume that radiation is the dominant excitation mechanism, this means that the 
outermost and star-facing portion of HD 189733b's atmosphere will be exposed to the most 
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radiation. This portion of the atmosphere will also be at the lowest density. Therefore, it is 
possible for some hydrogen to remain in an excited states at high altitudes, whereas at lower 
altitudes it will perhaps coUisionally de-excite. If the n = 2 hydrogen is only at the highest 
altitudes of HD 189733b's atmosphere, this may also explain its comparatively restricted 
velocity range — if it is already beyond the planet's Roche lobe (as we argued in §5.2.2p . it is 
not bound to the planet, and there is no physical or observational requirements that it be 
at a high velocity. The low density of this outer layer also means that it may be at a lower 
pressure than the rest of the atmosphere, consistent with the smaller velocity range. 



The best way to better understand the Ha abs orption in HD 18973 3 



be to explore its Lya absorption in more detail. IVidal-Madjar et al.l (|2003[ ) were able to 



j's ex osphere would 



study HD 209 458b's Lya absorption in great deta il due to their use of the E140M grating of 



STIS, whereas iLecavelier Des Etangs et al.l (120101 ) were unable to resolve the Lya absorption 



using ACS while STIS wa s not operational). Obsery ing HD 189733b's Lya absorption at 



a resolution comparable to IVidal-Madjar et al.l (120031 ) would enable a direct comparison of 



the velocity range over which each is present, and allow for a more detailed model of the 
Lya absorption (incorporating interstellar absorption, for example) and the resulting column 
density. 



5.3. The HD 209458b Feature 

5.3.1. Identification and Description 

In §4.31 we noted the unusual, strong feature in the transmission spectrum of HD 
209458b. We noted above that this feature is correlated with transit (or at least planetary 
phase) in our EMC results. All the reduction steps are performed "blindly" with no knowl- 
edge of whether an individual observation is in- or out-of-transit prior to the coaddition step; 
furthermore, we note that none of the other targets showed a similar feature in their master 
transit spectra. This argues that this feature is a real, planetary phase-correlated feature of 
unknown origin. However, we cannot completely rule out some unknown reduction effect, 
perhaps that is uniquely affected by some characteristic of the HD 209458b observations. 

We previously mentioned two broad categories of solutions that could be responsible 
for this feature — the Rossiter-McLaughlin Effect (RME), which creates a variation in the 
transmission spectrum that is dependent on wavelength, or a misalignment of the broad 
component of the Ha profile between the in- and out-of-transit spectra that is relative to 
the Ha line core and therefore not addressed by our process of higher-order wavelength 
correction. While these processes have the potential to match qualitatively what we see in 
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HD 209458b, we did not assess the quantitative plausibility of either scenario, which we do 
here. 



5.3.2. Rossiter-McLaughlin Effect? 

The RME is able to create features in a transmission spectrum that qualitatively match 
what we see — dips and/or spikes near stellar absorption lines that are the result of the planet 
blocking different portions of the star's spectrum in velocity spaced To test this, we made 
a simple simulation of the RME, in which we made a simulated Ha profile, assumed the 
stellar rotation speed, and calculated the RME based on the distribution of our observations 
during transit (i.e., which part of the velocity profile of the star is blocked in each observa- 
tion). Under a wide range of reasonable assumptions for this model (e.g., limb darkening 
parameters), we find multiple bases on which to reject the hypothesis that this large feature 
is the direct result of the RME. The first argument against this scenario is that the in-transit 
observations of HD 209458 are (randomly) weighted toward the first half of transit; when we 
simulate the RME, this results in spectral features with the opposite blue/red characteris- 
tics as the feature in the HD 209458b transmission spectrum. That is to say, features in the 
first half of transit block the blueshifted portion of the stellar disk, and should result in a 
"dip" to the blue and a "spike" to the red — the opposite of what we see. Furthermore, the 
feature we see in the HD 209458b transmission spectrum is too broad, by perhaps an order 
or magnitude or more, to be caused by the RME in a system where the star has a small 
known vsmi of only a few km s~^. Furthermore, in our simulations, the vertical magnitude 
of the feature (i.e., the height of the spike and dip) is also much smaller than the observed 
feature, again by a factor of several (the simulations would indicate that the magnitude of 
any RME is within the noise shown in Fig. |3]). Most importantly, we note that the RME is 
purely a transit effect, and we discussed in §4.31 and showed in Fig. [8] that the effect may be 
a more general phase effect with a maximum near transit. Thus we can rule out the RME 
as being responsible for this feature. 



5.3.3. Broad Stellar Ha Component Shifting? 

An alternative scenario for this feature is a relative velocity shift between the narrow 
and broad components of the stellar Ha line that is correlated with planetary phase. We 



^Note that this is superficially similar and related to, but nevertheless fundamentally different from, the 
feature that results from the RME in the more common plot of radial velocity as a function of time. 
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made a simple model of the Ha line by modeling it as two Gaussian components, one narrow 
and one broad. Using this method we fitted the master in-transit and out-of-transit spectra. 
When the narrow components of the master in- and out-of-transit spectra are aligned, the 
broad component of the out-of-transit spectrum is shifted by approximately —0.08 A (or 
—4 km s~^) relative to the in-transit spectrum. This model is a qualitative visual match to 
the spectrum, though the fits are statistically poor {xl ^ 1) due to the crudeness of this 
toy model. However, when we take the models obtained from these simple double- Gaussian 
model fits, one for the master in-transit spectrum and one for the master out-of-transit 
spectrum, and use them to create a synthetic transmission spectrum, we are able to replicate 
the observed transmission spectrum with x^ < 2. Thus it seems reasonably clear that this 
model can produce this feature in a way that, as described above in §5.3.2[ the RME cannot. 

We note that the —4 km s~^ shift is on the same order as the instrumental resolution of 
5 km s~ , though we are essentially performing a centroid measurement, which can be done 
with greater precision than the FWHM of an instrument's line spread function. This again 
raises the issue of a reduction artifact, which we discussed in §4.31 and §5.3. ![ but here we 
explore real, physical processes that might be responsible for such a shift that is strongly 
correlated with planetary transit. 

Some sort of SPl could conceivably cause this shift. As pointed out earlier, the evidence 
for SPl is very limited. We do note that highly interactive RS CVn binary systems display 
line behavior that is qualitatively similar to the feature in the transmiss ion spectrum of HP 
209458b (e.g., the observed vs. theoretical resid uals of Ha in IM Peg inlBiazzo et al.ll2006l ). 
Another extremely similar feature is shown in iHall &: Ramseyl (119921 ) in the "subtracted 
spectrum" shown for an eclipse of the AW Her system; this "subtracted spectrum" is the 
ratio of Ha to H/3. Neither of these cases a re perfect physi c al ana logues to a transit and 
our transmission spectrum — particularly the iHall &: Ramseyi ( 1l992l ) result, which is a ratio 
of two lines rather than a comparison of the same line in- vs. out-of-transit — but they do 
show characteristics that are similar to what we observe. 



The two primary ways in which a planet could influence its central star are through 
tidal and magnetic interactions. A tidal interaction could produce tidal bulges on the star , 



resulting in increased ste llar activity due to photospheric turbulence (JGuntz et al.l 12000 



Poppenhaeger et al.ll201ll ). Ma gnetic interac tions can increase stellar activity through in- 
creased magnetic reconnection ( jLanzal 120081 ) or h eating of the stellar atmosphere by flux 
tubes analogous to those in the Jupiter- lo system (ISchmittI 120091 ). If there is induced stellar 
activity on the surface due to the planet, it may produce emission that alters the shape of 
the broad or narrow component of the Ha as a function of whether the system is in a state o f 
in- vs. out-of-transit. For example, residual Ca II flux in HD 189733 (jShkolnik et al.ll2008[ ) 
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shows shifts in wavelength of up to ~0.25 A (see Fig. 4 of that paper). Perhaps an analogous 
residual flux in HD 209458's Ha line is of nearly constant intensity when averaged over the 
many observations but has a slightly shifted wavelength when the in-transit observations 
are compared to the out-of-transit observations. This shifting flux, imposed on the stellar 
line, might be able to create the narrow/broad relative shift that we see. Another possibility 
is that the magnetic influence of the planet directly influences the broad component of the 
stellar Ha line. If there is any physical separation in the layers of the stellar line — e.g., the 
core is dominated by a different layer than the wings — then this may be plausible. 

One difficulty with invoking stellar activity is that HD 209458 is not known to be 
particularly active, and is much less active than HD 189733. Why HD 209458, a less active 
star with a more distant planet (thus reducing the likelihood of SPI), would show this sort 
of effect is unclear. 



5.3.4. Shifting Due to an Existing Hydrogen Torus? 

A related alternative to the idea that SPI is causing the shifting of broad-component 
Ha in the star is the possibility that a component of the apparent stellar Ha line has a 
planetary origin. In other words, perhaps there a collection of material around the star, 
in a roughly toroidal shape, near the planet's orbit that formed from the escaping material 
from the planet. This would require that some of this hydrogen is excited and capable of 
exhibiting Ha absorption and that the planetary orbit affects the torus in some way related 
to its orbit. If this is a part of the broad component of the apparent stellar Ha profile, and 
it shifts in response to the planetary orbit, then it could be responsible for the feature that 
we see in the transmission spectrum. 

One difficulty with this interpretation is the extent of the torus. If it is only a partial 
torus near the planet, or a "comet-like" shape of material trailing the planet, then we would 
still expect to see an absorption feature in the master transmission spectrum. Given that 
we do not see such an absorption feature, for this general scenario of a torus of material 
to be correct would imply that the material exists continuously throughout the planetary 
orbit, and experiences a velocity shift in response to the planetary orbit. This would pose 
the problem of how the n = 2 hydrogen is created and sustained in the excited state. It 
would also be unclear why the shift would be similar at transit and at secondary eclipse, 
and different at quadrature. Further observations are needed to confirm or reject the feature 
that we see in HD 209458b, and, if confirmed, determine its nature and origin. 
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6. SUMMARY 

We perform atmospheric transmission spectroscopy of four exoplanetary targets in the 
spectral region surrounding the Ha hne. We find a strong feature of apparent Ha absorption 
in the transmission spectrum of HD 189733b. Our best attempts to remove stellar variability 
and assess other systematic errors indicate that this feature is real and correlated with 
planetary transit. The level of absorption is (—8.72 ± 1.48) x 10~^ integrated over a 16 A 
band and relative to the adjacent spectrum, significant at a 5.9cr level. Several complications 
exist in using this absorption to derive Te^c- We estimate Texc = 2.6 x 10^ after employing 
some simplifying assumptions, but discuss these assumptions and their effects in detail. In 
HD 209458b we find a dramatic feature with reflectional symmetry and discuss possible 
explanations. While we are unable to give a clear physical explanation for this feature, we 
discuss the clear correlation with planetary orbital phase. We place upper limits on the Ha 
absorption in our two remaining targets. These results all lend themselves to follow-up work 
to confirm and/or refine the properties of the observed features in HD 189733b and HD 
209458b, or to search for hydrogen absorption in other transitions of all four of these targets 
in order to better understand the nature hydrogen excitation in exoplanetary atmospheres, 
as well as possible insights into atmospheric escape in these planets. 
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-HD 156729 (telluric) 
-HD 147506 (target) 
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Fig. 1. — Examples of the spectra in the echelle order that contains H I A6562 for all four 
targets, along with the spectra of their respective telluric standard stars from the same night 
of observation. All spectra are reduced from the raw echelle images but are not processed 
with the additional steps we discuss in §2.21 (including removal of the Ha feature from the 
telluric spectra, telluric correction from the primary target, and normalization). The target 
spectra are shown on the count scale of the y-axis, but the telluric spectra are scaled to 
have the same maximum as the corresponding target spectra, and then arbitrarily shifted by 
25% of that same maximum. The global shapes of the spectra are dominated by the blaze 
function of the echelle. 
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Fig. 2. — A sample of a spline fit to one of the observations of the spectrum of star HD 
196867, the telluric standard star for HD 189733. The top panel shows the raw spectrum 
(scaled by 10^), along with its fit (thick red line), and the bottom panel shows the resulting 
residual, the normalized telluric absorption spectrum. 
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Fig. 3. — The transmission spectrum. The binned transmission spectrum is shown in red, 
with error bars (errors shown are the errors in the mean of the binned points). For reference, 
the master out-of-transit spectrum of the star is shown at the top green, scaled down to fit 
the plot, with a dotted green line showing the zero level and the top of each plot being unity 
in the normalized spectrum. A blue dashed line shows the zero point for the transmission 
spectrum. The purple vertical dot-dashed lines define the bandpass that is integrated to 
make our absorption measurements; from these lines to the edges of the plot define the 
comparison bandpasses. 
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Fig. 4. — Emprical Monte Carlot (EMC) analysis results. The distributions for the three 
versions of the EMC (see §3.2p are plotted in different colors, shown in the key: green is the 
"in-in" distribution, red is the "out-out" distribution, and blue is the "in-out" distribution. 
5000 iterations are used in each case. The master "in-out" measurement is shown as a 
vertical dot-dashed line. 
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Fig. 5. — The Lomb-Scargle periodograin of observations of HD 189733 Ha individual "trans- 
mission spectra" (i.e., individual out-of-transit spectra relative to the master out-of-transit). 
The horizontal red dashed-dotted line is the 0.3% (~ 3a) false alarm probability level. The 
dominant peaks are at ~ 1.086 and ~ 12.077 days, with many additional correlated peaks. 
The window function is also shown, with a peak at 0.9967 days, which is expected because 
the timing of the observations will create a peak at the sidereal period of 0.9972 days. 
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Fig. 6. — The individual out-of-transit transmission spectra of HD 189733 phased to the 
~ 12.077 day period, which corresponds closely to the stellar rotation period. The blue 
solid line is a low-order Fourier series fit to the out-of-transit points. The red dashed-dotted 
line is that same out-of-transit fit where we introduce an offset to fit the in-transit points. 
The systematic excess absorption seen in the in-transit data relative to the out-of-transit 
observations is evident. 
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Fig. 7. — Same as Fig. [6] but as a function of absolute time instead of phase. (No fit to the 
in-transit points is shown.) 
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Fig. 8. — This plot shows, for HD 209458, the individual transmission spectrum measure- 
ments. The Mabs here is for a 5 A band immediately to the blue of Ha line center, quasi- 
normalized to 5 A bands on either side. Therefore, a positive measurement is expected for 
what is seen in the master HD 209458b transmission spectrum, and a negative measure- 
ment for an inverted feature. Crosses are out-of-transit observations relative to the master 
out-of-transit spectrum; squares are in-transit observations relative to the same master out- 
of-transit spectrum. The top panel shows these measurements as a function of absolute 
time, while the bottom panel is phased to the orbital period of the planet. The color-coding 
of the top panel is kept for the bottom panel; note that the possible phase correlation in 
the bottom panel does not appear t o have any c o rrelat ion to date. The grey bar shows the 
phase range of the measurements by lWinn et al.l fl2004j ): the y- value of this bar is arbitrary. 
If the feature in the master transit spectrum has a physical origin (rather than an as-yet 
undetermined data reduction explanation), and if the feature is phase-correlated but not 
strictly transit- correlated (a s indic ated by the bottom panel), the limited phase range of the 
observations of IWinn et al.l (120041 ) compared to ours explains why we see this feature and 
they do not. 
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Fig. 9. — A plot showing how our derived T^xc varies as a function of n = 1 column den- 
sity A''i (based on the Lya absorption measurements), as well as (A^)LyQ/(A^)Ha. The 
top axis shows the 6- values (i.e., the Doppler parameter b = a\pl for a Gaussian velocity 
distribution) that correspond to these column densities. As 6 — t- oo, the absorption is unsat- 
urated and a lower limit on A'^i is reached; the left-hand axis is at 6 = 200 km s~ and nearly 
at this limit. Also shown are Tp^ (de fined in §5.2.5p and the limiting Tf,xc if the Lya W\ 



of iLecavelier Pes Etangs et al.l (120101 ) is an unsaturated measurement (which is unlikely). 
The red shaded area is the r ange of 6-values roughly implied by the fact that in Fig. 10 of 
Lecavelier Pes Etangs et al.l (120 lOl ) the ISM absorption is assumed to be saturated over the 
range ±25 km s~ . However, this is not a robust inference as a damped Lya line could show 
absorption at higher velocities even at low h. Note that if the Ha absorption is saturated, 
Tf^xc increases; if Lya absorption is underestimated due to the limits of the low-resolution 
ACS data (which is almost certain), then Tf,xc decreases. 



Table 1. Observations 



System 


In Transit'' 


Out of Transit^ 


Ingress/Egress^ 


Telluric 


Observations 




# of Obs. 


Time (s)'' 


# of Obs. 


Time (s)*^ 


# of Obs. 


Time (s)*^ 




# of Obs. 


<Time> (s)'' 


HD 147506'= 


28 


16508 


84 


49687 


5 


3000 


HD 156729'i 


34 


567 


HD 149026 


22(1) 


12718 (600) 


90 


52734 


5 


3000 


HD 156729"* 


34 


567 


HD 189733° 


38 (3) 


20140 (1640) 


173 (1) 


79011 (0) 


15(1) 


8243 (600) 


HD 196867 


45 


291 


HD 209458 


29(6) 


17310 (3510) 


72(3) 


43175 (1800) 


4 


2400 


HD 218045 


21 


80 



''Numbers shown are all observations taken, while the numbers in parentheses indicate observations that are discarded for various reasons, 
primary cloud-induced low S/N, but also various reduction artifacts. The "In Transit" and "Out of Transit" columns are mutually exclusive, 
defined by observation midpoint, while the "Ingress/Egress" observations are a subset of the "In Transit" observations, again defined by 
observation midpoint. 

''Times for primary science target exposures are total times, while telluric exposure times are average times, as the telluric observations are 
always used independently and are never coadded. 

"HD 147506 and HD 147506b are also known as HAT-P-2 and HAT-P-2b, respectively. 

"* Observations of HD 147506 and HD 149026 use the same telluric star. The values in the table this telluric represent the same observations, 
i.e. the 34 observations listed do not distinguish between which primary science target was originally associated with the telluric observation. 

"Some IRAF headers from the queue observing program prior to 12/2007 have "0 s" incorrectly listed in the exposure time header. This is 
true for 10 out of 211 observations of the HD 189733 system, and no observations for the other targets. Therefore, we estimate that the times 
shown for HD 189733 are incomplete by as much as 1800 s of in transit time (none discarded), 3000 s of out of transit time (600 s discarded, 
corresponding to the "0" in the table above), and 600 s of ingress/egress time (none discarded). 



O 



Table 2. System and Orbital Parameters of Targets 



Planet 



Reference Epoch 
(HJD-2400000) 



Period 

(days) 



Transit Duration 
(minutes) 



a/R. 



Rp/R. 



Inclination 

(°) 



References 



HD 147506b 
HD 149026b 
HD 189733b 
HD 209458b 



54387.4937 ± 0.00074 

54456.7876 ± 0.00014 

54279.43671 ±0.000015 

52854.825 ± 0.000135 



5.633473 ± 0.0000061 

2.87589 ±0.00001 

2.2185757 ± 0.00000015 

3.5247455 ± 0.00000018 



257.76 ± 1.87 
194.4 ±8.7 

106.1 ±0.7 

184.2 ±1.6 



;.99 



+0.39 
0.41 

+0.03 
0.81 



7.11] 

8.81 ±0.06 
8.76 ± 0.04 



0.07227 ± 0.00061 

U.U+yi_Q QQQ5 

0.15436 ±0.00022 
0.12086 ±0.0001 



90.0i°;» 
85.58 ±0.06 
86.71 ±0.05 



1, 2, 3, 4 
3, 5, 6 

3, 7, 8, 9 
3, 10 



Ref erenc es. — (1 _ 
|2003). ( 7)lAeol et al" 



Pal et al.l | |201(J). (2) iBakos et al.l 1I2OIOI'). (3) [Torres et al.l | |200^ 

' ' 2OI0I V (Sl lTriaud et al.l l l200g^ . (9) iBouchv et al.l l l2005fK (10) [Winn et al 






4')IWinn et al.l l l2007bl '). (5') ICarter et al.l ||2009| '). (6) IWolf et all 
2OO5I ). This tabic also made use of the Exoplanct Data 



Explorer jWright et al.ll201ll, |http://exoplanets.org|) and the Exoplanet Transit Database (http://var2.astro.cz/ETD/indcx.php). 
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Table 3. Ha Absorption Results 



Planet 



16 A Band 



6 A Band 



HD 147506b 
HD 149026b 
HD 189733b^ 



(-4.90 ±1.75) X 10"^ 

(3.07 ±1.53) X 10"^ 

(-10.6 ±2.43) X 10-^ 



HD 209458b (0.503 ± 0.802) x 10" 



-2.87 ±4.97) X 10"^ 
(7.76 ±4.20) X 10""^ 

-30.2 ± 6.33) X 10"^ 
(6.91 ±2.01) X 10"^ 



'^Uncalibrated values; see ^4.21 and Tabled 



Table 4. Additional HD 189733b Ha Absorption Results 



Velocity Range (km s 



Wx (ky 



Physical Significance 



-24.1,26.6] 
-59.5,59.5] 
-150,150] 
-365, 365] 



-143 (-118) 0.0155(0.0128) Range of Ho detection, tliis work 

-61.6 (-50.7) 0.0159 (0.0131) ii^esc at HD 189733b surface 

-27.5 (-22.6) 0.0180 (0.0148) Range (FWHM) of unresolved Lyo emission 

-10.7 (-8.8) 0.0171 (0.0141) Initial integration band of 16 A centered on Ha 



^Calibrated values in parentheses; see il4.2l 



